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We investigate the temperature region in which a Tomonaga-Luttinger liquid (TLL) de- 
scription of the charge sector of the one- dimensional Hubbard model is valid. By using the 
thermodynamic Bethe ansatz method, electron number is calculated at finite temperatures 
and fixed chemical potential. We observe maximum electron number as a function of tempera- 
ture close to the chemical potential of the upper critical value that corresponds to half filling. 
As the chemical potential approaches the upper critical value from below, the temperature 
(Tm) at which the electron number shows its maximum asymptotically approaches a universal 
relation. We show that, below the energy corresponding to Tm, the charge excitation spectrum 
nearly obeys a linear dispersion relation. The results demonstrate that Tm marks the important 
temperature below which TLL is realized. 

KEYWORDS: ID Hubbard model, electron number, excitation spectra, Tomonaga-Luttinger liquid, Bethe 
ansatz 



1. Introduction 

One-dimensional (ID) gapped many-body systems 
have attracted much attention. Interesting aspects have 
been revealed by the collaboration of theorists and ex- 
perimentalists. For ID gapped spin systems, it was ar- 
gued theoretically that the characteristic features of each 
system appear conspicuously in critical and dynamical 
properties above the critical magnetic field (^^c) where 
the gap is closed. Above the critical field, low-energy 
properties can be described as a Tomonaga-Luttinger 
liquid (TLL). The characteristic behavior depending on 
the model can be observed in, for example, the diver- 
gence of NMR relaxation rate with decreasing temper- 
ature.^' In fact, for the S = 1 Haldane-gap material 
(CH3)4NNi(N02)3^'''^ and the 5 = 1/2 bond-alternating 
chain material pentafluorophenyl nitronyl nitroxide^^ a 
divergence of NMR relaxation rate with decreasing tem- 
perature was observed. The field dependence of the di- 
vergence exponent was discussed in comparison with the- 
oretical results. ^"^^ Recent experimental advances have 
further allowed the observation of typical TLL behavior 
of thermodynamic quantities in ID gapped spin systems. 
The specific heat of the quasi- ID S = 1 bond-alternating 
antiferromagnet Ni(C9H24N4)(N02)C104 was measured 
at magnetic fields above Hc-^^ Low-temperature specific 
heat was found to be proportional to temperature, in- 
dicative of the gapless dispersion relation. The field de- 
pendence of its coefficient was in good agreement with 
the numerical result based on conformal field theory 
(CFT).9) 

Although the characteristic behavior of TLL has been 
studied intensively, it has been difficult to determine the 
temperature region where the TLL picture is valid from 
the gapped spin Hamiltonian itself. To develop precise 
comparison between the theoretical and experimental re- 
sults for the TLL. this issue has to be resolved. In a recent 



study, it was shown that the magnetization minimum ap- 
pears as a function of temperature close to H = He in 
ID gapped spin systems with axial symmetry. The 
magnetization minimum marks the important temper- 
ature below which TLL is valid. It was demonstrated 
further that the temperature for the magnetization min- 
imum approaches a universal relation as H approaches 
He from above. The arguments shed light on determin- 
ing the temperature region of the TLL by theoretical and 
experimental methods. Actually, the predicted features 
were confirmed theoretically for the S = 1 Haldane-gap 
system^°) and the ID S' = 1/2 Heisenberg model with a 
large Ising anisotropy^^^ in magnetic fields. 

The origin of the magnetization minimum and uni- 
versal relation was elucidated using the effective theory 
near the critical field and the low-temperature expansion 
of the free energy based on CFT.-'^"^ Since they are based 
on the fundamental nature of the TLL, the analyses seem 
to be applicable to general ID gapped systems with U(l) 
symmetry close to the critical point. To confirm this con- 
sideration, in this paper, we investigate the electron num- 
ber of the ID Hubbard model at finite temperatures us- 
ing the Bethe ansatz method. The charge sector of the 
ID Hubbard model is described by c = 1 CFT with 
U(l) symmetry. Here, c is the central charge. In 
§2, we outline the model and the method for the calcu- 
lation. In §3, we present the results. It is shown that a 
maximum appears in the temperature dependence of the 
electron number close to half filling. The temperature 
(Tm) for the maximum electron number approaches the 
universal relation, as the system approaches half filling. 
We find that, below the excitation energy correspond- 
ing to Tm, the linear dispersion relation of the charge 
excitation spectrum is nearly satisfied. Recently, vari- 
ous ID organic conductors have been synthesized. For 
such compounds, pressure- induced quantum phase tran- 
sitions from a charge-gapped state to the TLL state have 
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been investigated experimentally.^'*) The carrier number 
close to the critical point is expected to show a mini- 
mum/maximum as a function of temperature and the 
temperature for the minimum/maximum carrier number 
plays a role in determining the temperature region of the 
TLL. Section 4 provides a summary. 



2. Model and Method 

Let us consider the ID Hubbard model 

(i,j),c7 i i,<y 



(1) 



where the subscript (i, j) indicates the sum over nearest- 
neighbor sites, Cj^ annihilates an electron with spin 
at the iih lattice site, Ui^cr — and /i is 

the chemical potential. The hopping integral t is given 
in units of energy. As /i increases between its lower and 
upper critical values (/id = —2.0 < fi < ^ic2), electron 
number (n) increases monotonically from zero to unity. 
At half filling {n = 1), the Mott insulating state appears. 
The upper critical value of the chemical potential is ob- 
tained by the Bethe ansatz solution at zero temperature 
(T = 0) as /^c2 = F{Q)/D{Q)\q=^, where F{k) and D{k) 
are obtained using the integral equations 



2 cos k 



dk' cos fc'i?(sin k - sin k')F{k'),{2) 



D{k) = l+ dk' cos k'R{smk-smk')D{k'), (3) 

J-Q 

with k being the charge rapidity and R{x) = 
(l/27r) dwe"*'^^ (1 + e2'^l'^l)"\ The cutoff Q is ob- 
tained under the condition k{Q) = for the dressed en- 
ergy of charge rapidity. The dressed energy K{k) is ob- 
tained using the integral equation, 



K{k) 



2 cos k 



dXai{X — sinfc) 



dk' cos k'S{X- sink') K{k'), (4) 

where A is the spin rapidity, ai(A) ~ [U/'k)/{\^ + C/^), 
and S{x) = (l/4C/)sech(7rA/2f7). 

Thermodynamic quantities of the ID Hubbard model 
are formulated using the thermodynamic Bethe ansatz 
method. In the thermodynamic limit, an infinite set 
of nonlinearly coupled integral equations is derived for 
three types of pseudoenergies. Two of them are for two 
types of spin rapidities and the third one is for charge 
rapidity. The free-energy functional is expressed in terms 
of pseudoenergies. To perform the numerical calculation, 
we necessarily cut off two infinite sets of coupled integral 
equations at finite numbers, which are set at 60 and S.*®-* 
We employ the numerical technique developed by Usuki 
et aL,*^) where thermodynamic quantities are accurately 
calculated without numerical differentiation of the free- 
energy functional. The range of spin rapidity is set to be 
[—80,80], which is divided into 1600 points. The range 



of charge rapidity [— tt, tt] is divided into 800 points. The 
iterative calculations for pseudoenergies are carried out 
until they all converge within a relative error of 10^^. 
Their derivatives with respect to the chemical potential 
are then calculated in the same way within a relative 
error of 10~^. 

Using the thus-obtained pseudoenergies and their 
derivatives, we calculate electron number at finite tem- 
peratures. In the following calculations, we set U = 2.0, 
which leads to /ic2 = 1-66. 

3. Results and Discussion 



1 

0.98 
0.96 
0.94 
0.92 
0.9 





jU =1.700 — ^ 




1.625 — ^ 




1.600 




1.550 




, 1.500 









O.I 



0.2 



0.3 



0.4 0.5 



Fig. 1. (Color online) Temperature dependence of electron num- 
ber for a given chemical potential close to and above the upper 
critical value {fJ,c2 = 1-66). Electron number takes a maximum 
at the temperature indicated by an arrow. 
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Fig. 2. Temperature for maximum electron number as a function 
of chemical potential. The solid circles represent T^j and the 
solid line expresses the universal relation Tm = —xo{fi — /^c2), 
{xq ~ 0.76238). Inset: velocity as a function of chemical poten- 
tial, /ici = —2.0 and ij,c2 = 1-66. 



In Fig. 1, we show the temperature dependence of elec- 
tron number for several chemical potentials close to fj, — 
^c2- For the gapless regime < /j,c2, the electron num- 
ber shows a nontrivial maximum at low temperatures. 
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As chemical potential approaches the critical value fic2, 
Tm significantly decreases towards zero. For the Mott 
insulating region fi > /Xc2, by contrast, electron number 
shows no maximum and decreases monotonically with in- 
creasing temperature. From the low-temperature expan- 
sion for the free-energy functional based on CFT/^' 
electron number takes the form 

where hq is the electron number at T and Vc is the 
velocity of the charge excitation. The second term is the 
leading finite temperature correction. Accordingly, the 
behavior of the electron number near T = is deter- 
mined by the sign of dvc/dfj,. The velocity for the charge 
excitation is obtained as Vc = [dK{k)/dk]/[2'Kp{k)]\k=Q, 
where p{k) is the distribution function of the charge ra- 
pidity obtained using the integral equation, 

p(A) = — + cos/c / d\ai{X — smk) 



(a) 



2tt 



dk' S{X ~ sin k')p{k'). (6) 



As shown in the inset of Fig. 2, velocity decreases with 
increasing chemical potential close to p ~ ^,^2^ leading to 
a negative dvc/dp. As temperature is increased, there- 
fore, the electron number close to fi = fic2 increases near 
T — and then decreases, yielding a maximum structure. 

In Fig. 2, Tm is shown as a function of chemical 
potential. As chemical potential approaches /ic2, Tm 
asymptotically approaches the universal relation Tm — 
—xoin — A*c2), {xq ~ 0.76238). This relation can be de- 
rived similarly to that developed for ID gapped spin 
systems.^") Close to the upper critical value, the low- 
energy state can be approximated by E{p) ~ —p^ /2m + 
A*c2 ~ /^j where m is the effective mass of the elec- 
tron. The electron number for p < /ic2 is obtained 
using n = const. + ^mT/27rLii/2(-e(^'=^"'')/^) with 
Li„(a:) being the polylogarithm function. The differ- 
entiation of n with respect to T yields the condition 
for a maximum electron number, 2(/ic2 ~ p)/Tm — 
Lii/2(-e(''<^^-'')/'^")/Li_i/2(-e(^<^^~^)/^"), which leads 
to the universal relation mentioned above. If we replace 
p, with the magnetic field, the same universal relation as 
that derived for the ID gapped spin system^°) except for 
the sign can be obtained. The sign of the universal rela- 
tion depends on which critical point, the upper or lower, 
is paid attention to. The magnetization of ID spin sys- 
tems can be regarded as the particle number in fermionic 
particle language. Therefore, the thermodynamic quan- 
tity corresponding to the particle number shows a mini- 
mum/maximum as a function of the temperature close to 
the critical point of ID gapped many-body systems. The 
temperature for a minimum/maximum approaches the 
universal relation, as the system approaches the critical 
point. 

It is considered that the TLL is realized below Tm- 
We investigate the TLL behavior from the viewpoint 
of the excitation spectra. The charge excitation spec- 
tra of the ID Hubbard model are formulated using the 
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Fig. 3. (Color online) Lower edge of the charge excitation contin- 
uum for (a) /i = 1.48 and (b) ^ = 1.60 in the ID Hubbard model. 
The dotted lines are the excitation energies corresponding to Tm . 



Bethe ansatz solution. The excitation continuum can 
be obtained from the excitation energy and momentum: 
£(fch, fcp) = i^{kp) - K{kh) and p{kp, fch) = 2t: J^^^ p{k)dk, 
respectively, where —Q < kh < Q, — tt < kp < —Q, 
and Q < fcp < TT. In Fig. 3, we show the lower edge of 
the charge excitation continuum for two chemical poten- 
tials plotted in Fig. 2. The excitation energy (em) cor- 
responding to Tm for a given chemical potential is also 
shown in the figure by the dotted line. We find that for 
e < Em, the lower edge of the charge excitation contin- 
uum nearly obeys the linear dispersion relation, which 
is a characteristic of the TLL, although the deviation 
may be pronounced close to e ~ £m- As shown in Ap- 
pendix, also in the S — 1/2 Heisenberg model with a 
large Ising anisotropy in magnetic fields, the lower edge 
of the excitation continuum for H > He obeys the linear 
dispersion relation at e < e,„, where Em is the excita- 
tion energy corresponding to the temperature (Tm) for 
the magnetization minimum. In this way, the TLL for 
T < ThijTm in id gapped many-body systems has been 
confirmed from the spectral point of view. 



4. Summary 

We have investigated the temperature dependence of 
electron number in the ID Hubbard model using the 
Bethe ansatz method. We have found a maximum struc- 
ture at /i < /ic2. As chemical potential approaches /.tc2 
from below, Tm approaches the universal relation asymp- 
totically. We have also confirmed for the ID Hubbard 
model and ID S = 1/2 Heisenberg model with a large 
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Ising anisotropy that the hnear dispersion relation is 
satisfied below the excitation energy corresponding to 
Tim Tm- Judging from the findings obtained so far for ID 
gapped many-body systems, the thermodynamic quan- 
tity corresponding to the particle number shows a mini- 
mum/maximum as a function of the temperature close to 
the critical point and the TLL is realized at T < T„i, Tm- 
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Appendix: Excitation spectra of the S — 1/2 
Heisenberg model with a large Ising 
anisotropy in magnetic fields 
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Fig. A-1. (Color online) Lower edge of the excitation continuum 
for (a) H = 0.3 and (b) H = 0.4 in the S = 1/2 Heisenberg 
model with a large Ising anisotropy: A = 2.17. The dotted lines 
are the excitation energies corresponding to Tm. We use g/ig = 2. 
The antiferromagnetic exchange interaction is given in units of 
energy 



We calculate the excitation spectra of the S = 1/2 
Heisenberg model with a large Ising anisotropy in mag- 



netic fields using the Bethe ansatz solution. The Hamil- 
tonian is given by 

i 

i 

where J > and A — 2.17. We use gfiB = 2 and give J 
in units of energy. Excitation energy and momentum are 
expressed using the dressed energy 77(A) and the distri- 
bution function of the spin rapidity cr(A) as £(Ah, Ap) = 
77(Ap) — ?y(Ah) and p(Ap,Ah) — 27r J^^ a{X)dX, respec- 
tively, where —B < Xh < B, —n/(j) < Ap < —B, and 
B < Ap < 7r/(/) with (/) = cosh"^ A (0 > 0). The cutoff B 
is obtained from the condition ri{B) = 0, where r]{X) sat- 
isfies the integral equation 77(A) = 770(A) — jf^ dX'b2{X — 
X')ri{X') with 770(A) = 2H - (27rsinh0/0)6i(A) and 
6,1 (A) = sinh(7i0)/[cosh(770) — cos A]. The distribution 
function is obtained from the integral equation cr(A) = 
61(A) - J^g dX'b2{X - A')ct(A'). The lower edge of the 
excitation continuum and e,„ corresponding to Tm are 
shown in Fig. A-1. We also find that in this model, the 
lower edge of the excitation continuum obeys the linear 
dispersion relation at e < Em. 

1) R. Chitra and T. Giamarchi: Phys. Rev. B 55 (1997) 5816. 

2) N. Haga and S. Suga: J. Phys. Soc. Jpn. 69 (2000) 2431. 

3) T. Goto, Y. Fujii, Y. Shimaoka, T. Maekawa, and J. Arai: 
Physica B 284-288 (2000) 1611. 

4) T. Goto, T. Ishikawa, Y. Shimaoka, and Y. Fujii: Phys. Rev. 
B 73 (2006) 214406. 

5) K. Izumi, T. Goto, Y. Hosokoshi, and J. -P. Boucher: Physica 
B 329-333 (2003) 1191. 

6) T. Sakai and M. Takahashi: J. Phys. Soc. Jpn. 60 (1991) 3615; 
Phys. Rev. B 43 (1991) 13383. 

7) S. Suga and T. Suzuki: Physica B 346-347 (2004) 55. 

8) T. Suzuki and S. Suga: Phys. Rev. B 70 (2004) 054419. 

9) M. Hagiwara, L. P. Regnault, A. Zheludev, A. Stunault, N. 
Metoki, T. Suzuki, S. Suga, K. Kakurai, Y. Koike, P. Vorder- 
wisch, and J. H. Chung: Phys. Rev. Lett. 94 (2005) 177202. 

10) Y. Maeda, C. Hotta, and M. Oshikawa: Phys. Rev. Lett. 99 
(2007) 057205. 

11) S. Suga: J. Phys. Soc. Jpn. 77 (2008) 074717. 

12) N. Kawakami and S.-K. Yang: Phys. Lett. A 148 (1990) 359. 

13) H. Frahm and V. E. Korepin: Phys. Rev. B 42 (1990) 10553. 

14) See, for example, L. Degiorgi and D. Jerome: J. Phys. Soc. 
Jpn. 75 (2006) 051004. 

15) M. Takahashi: Prog. Theor. Phys. 47 (1972) 69. 

16) There is an alternative method for calculating the exact ther- 
modynamics of the Hubbard model, G. Jiittner, A. Kliimper, 
and J. Suzuki: Nucl. Phys. B 522 (1998) 47169. In this method, 
an infinite set of nonlinearly coupled integral equations is 
avoided fundamentally and no truncation is performed. 

17) T. Usuki, N. Kawakami, and A. Okiji: J. Phys. Soc. Jpn. 59 
(1990) 1357. 

18) H. W. J. Blote, J. L. Cardy, and M. P. Nightingale: Phys. Rev. 
Lett. 56 (1986) 742. 

19) I. Affleck: Phys. Rev. Lett. 56 (1986) 746. 

20) F. H. L. Essler, H. Frahm, F. Gohmann, A. Kliimper, and V. E. 
Korepin: The One- Dimensional Hubbard Model (Cambridge 
University Press, Cambridge, 2005) pp.209-267. 



